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We have studied the trends in CO adsorption on close-packed metal surfaces: Co, Ni, Cu from the 
3d row, Ru, Rh, Pd, Ag from the 4d row and Ir, Pt, Au from the 5d row using density functional 
theory. In particular, we were concerned with the trends in the adsorption energy, the geome- 
try, the vibrational properties and other parameters derived from the electronic structure of the 
substrate. The influence of specific changes in our setup such as choice of the exchange correla- 
tion functional, the choice of pseudopotential and size of the basis set, substrate relaxation has 
been carefully evaluated. We found that while the geometrical and vibrational properties of the 
adsorbate-substrate complex are calculated with high accuracy, the adsorption energies calculated 
with the gradient-corrected Perdew-Wang exchange-correlation energies are overestimated. In ad- 
dition, the calculations tend to favour adsorption sites with higher coordination, resulting in the 
prediction of wrong adsorption sites for the Rh, Pt and Cu surfaces (hollow instead of top). The 
revised Perdew-Burke-Erzernhof functional (RPBE) leads to lower (i.e. more realistic) adsorption 
energies for transition metals, but to wrong results for noble metals - for Ag and Au endothermic 
adsorption is predicted. The site preference remains the same. We discuss trends in relation to the 
electronic structure of the substrate across the Periodic Table, summarizing the state-of-the-art of 
CO adsorption on close-packed metal surfaces. 

PACS numbers: 



I. INTRODUCTION 



The development of modern theoretical surface science 
provides an opportunity to investigate surfaces and ad- 
sorbate structures on the atomic scale with useful appli- 
cations in industrial technologies. Much effort has been 
devoted to study CO chemisorption and dissociation on 
transition metals. There arc numerous papers and re- 
views which deal with this system from different points 
of view (electronic, structural, vibrational) 1, 3,11, 0|- 
One of the central questions is how the strength of the 
chemisorption of CO and the preference for the specific 
adsorption site varies across the transition metal (TM) 
series. A particular case, CO adsorption on the Pt(lll) 
surface, has attracted much attention in the past, since 
for this system state-of-the-art DFT calculations fail in 
predicting the correct site preference 0. The question 
which immediately arises, whether this case is an excep- 
tion or the rule ? 

Additionally, a vast number of theoretical papers ap- 
peared in the literature since Ying, Smith and Kohn 
presented in the mid 70's a first self-consistent density 
functional study of chemisorption on metal surfaces (H 
on tungsten) [fl] ■ The number of theoretical adsorption 
studies of TM molecules on different surfaces is increas- 
ing with time not only because of their importance in 
catalysis, but also due to the increasing reliability of the 
"measured" properties. In the past, several systematic 
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studies of the adsorption of CO molecules on transition 
metal surfaces have been performed I'^ 1990 

N0rskov proposed a model of chemisorption on transi- 
tion metal surfaces which was later expanded and is now 
quite generally accepted [3> 0- A main feature of the 
model is the importance of the position of the d-band 
center relative to the HOMO and LUMO of the adsor- 
bate. The importance of understanding the correlation 
between the geometric and the electronic structure arises 
from the proposed mechanistic model for chemisorption. 
One among several trends is the correlation between the 
CO stretching frequency vc-o and the energy level dif- 
ference between ha and Itt orbitals of the adsorbed CO 
(A(5(T — Itt)) proposed by Ishi et al. jlflj. 

In this paper we present an extensive density func- 
tional study of the adsorption of CO on the closed-packed 
surfaces [(111) for face-centered cubic, resp. (0001) for 
hexagonal metals] of Co, Ni, Cu, Ru, Rh, Pd, Ag, Ir, Pt 
and Au. After a description of the set-up in Section ITU 
we characterize briefly the clean surfaces in Section mil 
Section im is devoted to the adsorption of CO. Starting 
from the geometric structure of the adsorbate-substrate 
system, going through vibrational and electronic prop- 
erties (stretching frequencies of the adsorbate, and of 
the adsorbate-substrate bond, occupation of the anti- 
bonding 27r*-like orbital, density of states, redistribution 
of the charge density) as well as site preference we draw 
a complete picture of the CO adsorption on the close- 
packed TM surfaces. Moreover, we investigate the influ- 
ence of the exchange-correlation functional and the cut- 
off energy on the site preference (Section 0). In Section 
fVllwe discuss our results in the light of the experimental 
literature and analyze trends and correlations between 
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the investigated properties. This paper tries to go beyond 
a pure table of references from experiments and different 
DFT calculations, by providing a large and consistent 
database, in which each element was treated in exactly 
the same manner. Further, one of our main goals is not 
only to obtain theoretical values of spectroscopic accu- 
racy, but also to derive useful trends of CO adsorption 
on 3d, 4d and 5d transition metal surfaces with the hope 
of applicability in the prediction of the adsorption and 
catalytic behavior. 



II. METHODOLOGY 

The calculations in this work are performed u sing 
the Vienna ab-initio simulation package VASP 0, Il2l | 
which is a DFT code, working in a plane-wave basis 
set. The electron-ion interaction is described using the 
projector-augmented- wave (PAW) method [H [l| with 
plane waves up to an energy of Ecut = 450 eV (for some 
calculations harder pseudopotentials were used for C and 
O which require energy cut-off of 700 eV) . For exchange 
and correlation the functional proposed by Perdew and 
Zunger 15] is used, adding (semi-local) generalized gra- 
dient corrections (GGA) of various flavor (PW91 [iq . 
RPBE 01). These GGAs represent a great improve- 
ment over the local density approximation (LDA) in the 
description of the adsorption process. 



sampled by a grid of (4x3x1) k-points. We have chosen 
this coverage as a compromise between small adsorbatc- 
adsorbate interactions ("low coverage limit") and low 
computational effort. 

In the calculation we investigated the adsorption on 
the close packed surfaces of 10 metallic elements from the 
3d (Co, Ni, Cu), 4d (Ru, Rh, Pd, Ag) and 5d transition 
metal rows (Ir, Pt, Au) of the Periodic Table. For all 
hcp-elements, Co and Ru, the (0001) surface was used 
with the ideal c/a ratio of 1.63. The spin-polarization of 
Ni and Co was also taken into account. 

Vibrational properties of CO were computed by ap- 
plying a finite-differences method to create the Hessian 
matrix which we diagonalize to obtain the characteristic 
frequencies. We have calculated the metal-CO (j^m-co) 
and C-0 stretch (i^c-o) frequencies in the direction per- 
pendicular to the surface plane. 

The free CO molecule is characterized by a calculated 
stretch frequency of 2136 cm~^ at an equilibrium bond 
length of 1.142 A. The corresponding experimental val- 
ues are 2145 cm~ 1 and 1.128 A^. The problem of 
a too large CO binding energy {E^q^^^ — 11.76 eV, 
^co ~ 11.45 eV |0|) stems mainly from the error in 
the energy of the free atom, where high density gradients 
make an accurate description more difficult. 

Further, the bonding of adsorbate to the surface by 
calculating density of states and charge density flow was 
investigated. 
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FIG. 1: Top and side view of the slab used in the calculations. 
In the top view of the c(2x4) cell we denote all investigated 
high-symmetry sites. 

The substrate is modelled by four layers of metal sep- 
arated by a vacuum layer of approximately double thick- 
ness, as shown in Fig. ^ The two uppermost substrate 
layers and the CO molecule are allowed to relax. This 
enables us to check the influence of the relaxation on the 
adsorption system. The Brillouin zone of the c(2x4) sur- 
face cell (equivalent to a coverage of 8 = 0.25 ML) was 



III. BULK AND CLEAN SURFACES 
A. Lattice constant 

For a general understanding of the adsorption process, 
especially for the adsorption in higher coordinated sites, 
the lattice constant of the substrate is an important pa- 
rameter. The optimal adsorption height for example will 
always be determined by an interplay between the op- 
timal carbon-metal bond-length and the lattice param- 
eter. For that reason, we give in Fig. [3 the theoretical 
lattice constants together with the experimental values 
[20l |. These lattice parameters are compared with lat- 
tice parameters calculated with the PW91 and RPBE 
exchange correlation functional. 

The overestimation of the lattice parameter (calculated 
within the GGA) is characteristic for the heavier ele- 
ments. The difference between the experimental and cal- 
culated lattice constant is around 1% for 3-d row metals 
and around 2% for other elements. Such a 2% difference 
in the lattice constant corresponds to a change in the CO 
adsorption energy of about 0.03 eV for Pd(lOO) 21] and 
Ru(OOOl) surface jl^. In general, the lattice parame- 
ter increases along the rows and columns of the Periodic 
Table. For a greater lattice constant one might expect 
that the adsorbate will come closer to the surface. On 
the other hand, as the width of the d-band increases, the 
binding is reduced and the actual height of the adsorbate 
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FIG. 2: Lattice parameters for a part of the Periodic Table. 
We show the experimental and calculated lattice parameters 
for PW91 and RPBE functional. Although, Ru and Co are 
hep metals we have included them together with fee metals 
at an ideal c/a = 1.63. 



over the surface should increase. 



Work-function 



In Fig. 01 we present the calculated work- functions to- 
gether with experimental values. 
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FIG. 3: Experimental and calculated values of the work- 
function (PW91) for various transition metal (111) surfaces. 
References for experiments: Co |23I . Ni 23], Cu |23, Ru [2^ . 
Rh Pd Ag Ir |3, Pt'lMlil, Au 

We find that the calculated work- functions (PW91) are 
in all cases slightly lower (by at most 6% for Ag) than the 
measured values. The largest discrepancies in terms of 
absolute values are found for 4d metals, where the differ- 
ences for Ru, Rh and Ag surface are wO.3 eV. However, in 
general the agreement is quite good. The work-function 



always increases with the d-band filling. Only for the 
noble metals (Cu, Ag, Au) it decreases again. Similarly, 
the work-function increases along the columns when go- 
ing down from the 3d to the 5d metals. 



C. Electronic structure : d-band center 

The position of the d-band center of the clean surface is 
another important characteristic which is closely related 
to the strength of the CO-surface interaction. As it was 
already argued 0, [s^l , the d-band centers play a sig- 
nificant role in the bonding for many adsorbate-substrate 
systems where the major interaction is due to hybridiza- 
tion of the HOMO and LUMO of the adsorbate and the 
d-orbitals of the substrate. 
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FIG. 4: The d-band centers of the clean TM surfaces inte- 
grated until Epcrmi Icvel and 15 eV above Epcrmi, where at 
least 9.5 electrons are accommodated in the d-band. As a 
comparison the calculated d-band centers from the article of 
Hammer and N0rskov are taken lal. 



However, simple as it might seen to be, the proper def- 
inition of the d-band center is not that straightforward: 

(i) Within a plane-wave method, partial angular- 
momentum decomposed densities of states can be defined 
in different ways. We chose the projection of the plane- 
wave components onto spherical orbitals within some 
atomic radius. 

(u) DFT apphes only to the ground-state, excited 
states are usually predicted at too low energies. There- 
fore we report in Fig. 01 the center of gravity of the occu- 
pied d-band (integration up to the Fermi- level) d^, and 
the center of the entire d-band d;!.^ (integration from the 
lower band edge up to 15 eV above EFcrmi where at least 
9.5 electrons have been accommodated in the d-band), 
together with the d-band centers taken from the older 
LMTO calculations of Hammer and N0rskov (df ^) 0. 
The d^ and d;!.^ values provide upper and lower bounds 
for our d-band centers. The d^^ lies in between for TM 
surfaces, but it is lower in energy for noble metals. This 
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stems from a small contributions to the density of states 
far from the main d-band peak. As per definition, the 
difference between d° and dj^ increases with decreasing 
band-filling, d^^ shows only little variation in groups 
VIII to IB. As expected, it drops significantly for the 
noble metals compared to the metals in the same row. 



IV. CO ADSORPTION 

We have studied the adsorption of CO in the top, 
bridge and hollow (fee, hep) sites in the c(2 x 4) cell, 
corresponding to a coverage of a quarter of a monolayer 
(6=0.25 ML) (see Fig. From these calculations we 
derive in the following trends in the adsorption energies, 
the vibrational properties, the geometric and the elec- 
tronic structure. 



A. Adsorption energy and site preference 

Beside the adsorption geometry and the vibrational 
properties a precise determination of the adsorption en- 
ergy for a given geometry is one of the keys to under- 
standing the mechanism of the catalytical activity and 
selectivity [sj- The adsorption energy describes the 
strength of the chemical bond between the CO molecule 
and the metallic surface. 

The adsorption energy for the stable adsorption site 
corresponds to the absolute minimum on the potential 
energy surface (PES) of the molecule moving on the sur- 
face. The adsorption energies calculated for other high- 
symmetry sites correspond either to local extreme or to 
a saddle point on the PES. The position and the height 
of the saddle points define the diffusion path and the 
activation energy for surface diffusion. 

The experimental isosteric heat of the adsorption 
Hads(0) is not directly comparable with the calculated 
adsorption energy. If we want to compare our theoret- 
ically calculated adsorption energy Eads (i-e. the gain 
in energy when a certain amount of CO is adsorbed) at 
0=0.25 ML with the isosteric heat of adsorption (a dif- 
ferential energy gained when adsorbing one additional 
molecule on a surface at a certain CO coverage) we need 
to take the integral up to the coverage used in the calcu- 
lation: 



^?ad.,int(ecalc) = / ilads (e)de. (1) 

f calc Jo 

The experimental heats of adsorption with and the pro- 
posed adsorption sites are presented in TableQ] We show 
adsorption sites which are found at 0.25 ML or up to 
0.5 ML coverage. 

The main trends in the experimental results given in 
Table m include: 



1. CO adsorption on TM surfaces is several times 
stronger (Eads ~ 1-2 to 1.7 eV) than on noble met- 
als (Eads 0.3 to 0.5 eV). 

2. All measured adsorption sites (with the exception 
of on-top site on Ru(OOOl)) show a slight decrease of 
the heat of adsorption going from the zero coverage 
limit to a quarter-monolayer coverage. 

3. The heat of adsorption decreases with increasing 
filling of the d-band. 

Fig. [51 summarizes the calculated adsorption energies 
(using the PW91 GGA) for the high symmetry positions 
described in Fig. ^ For the transition metals we note 
a pronounced tendency to overestimate the adsorption 
energies - this is independent of an eventual relaxation 
of the surface. For the noble metals (with the exception 
of Cu) on the other hand, the calculations rather tend 
to underestimate the adsorption energies - this could be 
related to the neglect of dispersion forces. 

Experimentally, the adsorption energies on the TM 
surfaces are found to decrease monotonously with in- 
creasing band filling. In the calculation, the variation 
depends on the adsorption site: Eads decreases for on-top 
adsorption, but increases for sites with higher coordina- 
tion. This leads to wrong predictions for Cu, Rh and 
Pt where the hollow sites were found to be preferred for 
CO adsorption, whereas experiment shows that adsorp- 
tion occurs on the top sites. Moreover, the relaxation of 
the metallic substrate is necessary for Co and Ru sur- 
faces, i.e. if the surface is not allowed to relax we get 
a different site prediction compared to the experiments. 
However, even after substrate relaxation, the energies for 
CO adsorption in hollow and top sites on Co, Ag, and 
Au surfaces remain almost degenerate. Site preference 
is indeed a fundamental problem for CO adsorption on 
metallic surfaces. We shall try to elucidate the origin of 
this problem below. 

Another interesting point is that among the hollow 
sites, the preference changes from hep (Co, Ni, Ru, Rh, 
Ir) to fee (Cu, Pd, Ag, Pt, Au) with increasing d-band 
filling. 

B. Geometric structure 

In the past many studies have demonstrated that on 
the basis of DFT a reliable determination of the geo- 
metrical structure is possible. The local adsorption ge- 
ometry of the CO molecule on an unrelaxed TM sur- 
face can uniquely be characterized by the S(urface)-C 
distance (height above the surface) and the C-0 bond 
length. If additionally the metal atoms are allowed to re- 
lax, this modification of the surface can be described by 
the change of the inter-layer distance and the buckling of 
the surface. 

On all studied surfaces the CO molecule adsorbs in an 
upright position with the carbon atom pointing towards 



5 



Experiment 


Theory 


surface 


site 


Hads(OML) H 


ads (i ML) H 


adSjint 


method 


e [ML] ref. 


■p\j-'wyi,45u 

ads 


ads 


ads 


Co 


top 


-1.33 
-1.19 


-1.33 


-1.33 


SP 
TDS 


0.33 ^ 

[33] 


-1.65 


-1.62 


-1.32 


Ni hollow 


-1.35 


-1.28 


-1.32 


CAL 


0.25 [34] 


-1.95 


-1.90 


-1.44 






-1.31 


-1.28 


-1.30 


TDS 


0.25 [35] 












-1.31 


-1.06 




TDS 


0.25 [36] 








Cu 


top 


-0.49 
-0.425 
-0.46 






TDS 
TDS 
TDS 


[37] 
[38] 
0.33 [39] 


-0.75 


-0.72 


-0.42 


Ru 


top 


-1.66 


-1.66 


-1.66 


TDS 


0.20 [40] 


-1.89 


-1.82 


-1.69 


Rh 


top 


-1.5 


-1.38 S: 


-1.44 


TDS 0.2-0.25 [411 


-1.89 


-1.86 


-1.55 






-1.37 


-0.77 a 


-1.07 


TDS 


0.25 [42j 












-1.71 






TREELS 


[43] 












-1.40 


-1.28 


-1.36 


He 


[44] 








Pd hollow 


-1.30 






TDS 


[45] 


-2.14 


-2.09 


-1.68 






-1.54 


-1.30 a 


-1.42 


TDS 


0.25 [46] 








Ag 


top 


-0.28 


TDS 


[47] 


-0.16 


-0.14 


0.18 


Ir 


top 


-1.81 


-1.55 
-1.52 


-1.63 


TDS 


0.25 [48] 
0.33 [49] 


-1.98 


-1.94 


-1.64 


Pt 


top 


-1.50 
-1.39 
-1.43 


-1.22 
-1.30 ft 


-1.33 
-1.37 


TDS 
LITD 
TDS 


0.25 [50] 
0.25 [51] 
0.25 [52] 


-1.70 


-1.67 


-1.34 


Au 


top 


-0.40 


TDS 


[53] 


-0.32 


-0.24 


0.12 



TABLE I: Experimental heats of adsorption for 0=0 and 0.25 ML (if available) and the integral heat of adsorption between 
zero limit and 0.25 ML. Integral heats of adsorption Hads.int are compared with our calculated adsorption energies Eads at the 
experimentally observed adsorption sites for the PW91 and RPBE exchange-correlation functionals using two different sets 
of pseudopotential (Ecut = 450 eV, Ecut ~ 700 eV, cf. text). TDS - Thermal Desorption Spectroscopy, CAL - Calorimetry 
method, He - Helium scattering method, LITD - Light Induced Thermal Desorption, SP - SP measurement, TREELS - Time 
Resolved Electron Energy Loss Spectroscopy. 
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FIG. 5: Calculated CO adsorption energies on TM surfaces 
for the PW91 exchange-correlation functional and energy cut- 
off 450 eV. Experimental CO heat of desorption are labelled 
by the stars (see Table 0. 



FIG. 6: Calculated heights of the CO adsorbate (surface-C 
distance) in different adsorption sites on various TM surfaces 
together with values taken from experimental literature (Ta- 
ble EJ. 



the surface. The tilt with respect to the surface normal 
is always less than 2°. 

The minimum distance between neighboring adsor- 
bates varies in our setup (c(2x4)) between 4.3 A (Co) 
to 5.1 A (Au). Therefore direct as well as indirect inter- 



actions can be considered to be small and our calculations 
probe the low-coverage case. 

In response to the adsorption of the CO molecule, the 
metal atoms in the immediate surrounding move out- 
wards. This effect is quite localized. The buckling is 
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less pronounced for adsorption in higher coordinated sites 
(~ 0.1 A) and increases for one-fold coordinated CO in 
an on-top position to ^0.2 A. 

The height of the CO molecule above the surface 
(ds-c) is determined by the extension of the metal d- 
orbitals and (for the higher coordinated sites) by the 
nearest neighbor distance of the substrate atoms. The 
greater the extension of the interacting metal orbitals the 
larger is the distance between the C-atom and the closest 
metal atom (dM.^c)- Similarly, for a given metal-C dis- 
tance, the molecule will adsorb at lower height for greater 
lattice constants. The essential importance of the metal- 
carbon bond distances for understanding the mechanism 
of the metal-CO bonding was pointed out earlier ^3 ■ 




Co Ni Cu Ru Rh Pd Ag Ir Pt Au 



FIG. 7: Calculated C-0 bond lengths in top, bridge, fee and 
hep hollow sites on TM surfaees together with experimental 
values as listed in Table HTl 

As it was already indicated in the previous studies, 
a different occupation of the interacting d-band orbitals 
leads to different values for the ds-c and dM-c lengths. 
In Figs. El and [3 the main characteristics of the CO ad- 
sorption geometry on TM surfaces are compiled. 

When comparing the adsorption geometry, the noble 
metals behave differently from all other metals. This re- 
flects the completely filled d-band and hence the different 
bonding mechanism. Adsorption on the noble metals is 
significantly weaker (see adsorption energies Fig. [SJ, and 
consequently the ds-c and dM-c lengths are enhanced 
(see Fig. 0). On the other hand, CO adsorption on the 
surfaces of TM is stronger and ds_c (resp. dM-c) tends 
to decrease as the filling of the d-band increases. This 
trend is more pronounced for 4d than for 5d or 3d met- 
als (for Co and Ni this does not hold for on-top adsorp- 
tion). For on-top adsorption, the height of the adsorbate 
increases going from 3d via 4d to 5d elements. The co- 
ordination has a strong effect on the CO height: CO in 
the lower coordinated sites is closer to the surface. 

The C-0 bond length (dc-o) of the free molecule 
was calculated to be 1.142 A. On adsorption, the dc-o 
length increases with the coordination of the adsorp- 



tion site, similarly to ds-c and dM-c: dc-o rises from 
1.15 - 1.17 A for one-fold to 1.17 - 1.20 A for three-fold 
coordinates sites. This trend was already stressed in the 
experimental work of Westerlund et al. [s^ . Addition- 
ally, the C-0 bond is more elongated as the filling of the 
d-band decreases. 

There are many data for the geometrical structure of 
the CO adsorption on TM from past experiments 55]- 
[6^ . Only for iridium, gold and silver surfaces we are 
not aware of any experimental literature. A valuable col- 
lection of the detailed experimental data for the adsorp- 
tion of small molecules on the metallic surfaces can be 
found in the study of Over Q . A part of the vast data is 
presented in Table ITU 

The main difficulty in comparison with experiment 
(LEED, XPD, EXAFS, ...) is that the estimated error for 
the C-0 bond length (~ ±0.05A) is mostly too large, cov- 
ering the whole interval of all calculated dc-o values. Al- 
though the c(2x4)-CO structure used in the calculation is 
not realized on all metallic surfaces, values for p(\/3 x \/3) 
and p(2 x 2) structures give values in line with our re- 
sults. This implies that the CO adsorbate-adsorbate in- 
teractions at a coverage lower than O = ^ ML play only 
a minor role. 



C. Vibrational frequencies 

Vibrational frequencies can be measured very accu- 
rately using infrared spectroscopy (RAIRS) , electron en- 
ergy loss spectroscopy (EELS) and sum frequency gener- 
ation (SFG) methods. The dependence of the stretching 
frequency vc-O on the coordination by the surface atoms 
has often been used as an indication of the adsorption 
site. The correlation between the vibrational and elec- 
tronic properties was analyzed in the paper by Ishi et al. 
llOj l in which a correlation between I'c-o and the energy 
level difference between So- (5cr character after adsorp- 
tion) and Itt orbitals of the adsorbed CO (A(5i7 — Itt)) 
was suggested. 

In the following we will focus on two eigenmodes: the 
M-CO (t'M-co) and the C-0 {vc-o) stretching frequen- 
cies. Our calculated I'c-o values are compiled in Fig. 
IHl For free CO molecules a stretching frequency of vc-O 
2136 cm"^ was calculated, this value is lowered due to 
the formation of the bond between the molecule and the 
metallic surface. The C-0 vibrations exhibit the ex- 
pected strong coordination dependence: CO molecules 
adsorbed in the lower coordinated sites vibrate faster, 
typical values for the C-0 stretching frequency are 1990- 
2100 cm-i, 1830-1880 cm-i, 1750-1810 cm^^ for carbon 
monoxide in top, bridge and hollow sites. The stretching 
frequency for the CO molecule in bridge and hollow sites 
on noble metals is higher than on the TMs whereas for 
on-top adsorption the frequencies are almost the same 
for the noble metals and those of the Pt group. 

The dependence of vq-o on the d-band filling of the 
substrate is also obvious: increased d-band filling raises 
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Experiment 


Theory 


surface 


site ds-c (A) dc-o (A) 


e (ML) ref. 


d^s^J^c (A) dg!! 


o (A) 


Co 


top 1.78±0.06 1.17±0.06 


0.33 [55] 


1.83 


1.166 


Ni 


hollow 1.34±0.07 1.15±0.07 


0.5 [56] 


1.329 


1.192 




hollow 1.29±0.08 1.18±0.07 


0.5 [56] 


1.327 


1.192 


Cu 


top 1.91±0.01 


0.33,0.44 [57] 


1.96 


1.156 


Ru 


top 1.93±0.04 1.10±0.05 


0.33 [58] 


2.03 


1.166 


Rh 


top 1.87±0.04 1.20±0.05 


0.33 [59] 


1.99 


1.162 


Pd hollow 1.27±0.04 1.14+" !,' 


0.33 [60] 


1.31 


1.188 




hollow 1.29±0.05 1.15±0.04 


0.33 [61] 


1.33 


1.188 


Pt 


top 1.85±0.10 1.15±0.1 


0.3 [62] 


2.00 


1.157 



TABLE IL Experimental geometrical structures from low energy electron diffraction (LEED) studies of the CO adsorption on 
TM surfaces together with the values obtained in this study (superscript "cal"). The normal distance between carbon and the 
TM surface (ds-c) and C-O bond length (dc-o) at the coverage Q together with corresponding references. 




Co Ni Cu Ru Rh Pd Ag Ir Pt Au Co Ni Cu Ru Rh Pd Ag Ir Pt Au 



FIG. 8: Calculated C-O stretching frequencies (i/c-o) on 
closed-packed TM surfaces together with experimental values 
as listed in Table IIIII The calculated value of vc-o for the 
free CO molecule is 2136 cm~^ compared to the experimental 
value of 2145 cm~^. 



i/C-o; the correlation is essentially linear. An exception 
to the rule is CO in on-top sites on noble metals. With 
decreasing d-band filling the stretching frequencies for 
the CO adsorbed in the fee and hep hollows begin to 
differ, with the lower vc-O for the hep sites. 

The adsorbate-substrate stretching frequencies vm-co 
are compiled in Fig. |3 j^m-co exhibits a similar coordi- 
nation dependence as i^c-o- CO molecules adsorbed in 
the lower coordinated sites vibrate faster, typical inter- 
vals of the vm~co stretching frequencies on the TM sur- 
faces are 400-475 cm~\ 340-380 cm~\ 300-350 cm^^ for 
top, bridge and hollow sites. There is no pronounced cor- 
relation with the d-band filling. For noble metals vm-co 
lies generally below 300 cm^^ and there is only a very 
weak site dependence. Once more, this reflects the weak 
CO-noble metal bonding. The metal-substrate stretching 
frequency is particularly low for the Ag surface, consis- 
tent with a very low adsorption energy. We have col- 
lected available experimental data for vc-o smd vm-co 



FIG. 9: Calculated M-CO vibrational frequencies (vm-co) 
on closed packed TM surfaces together with experimental val- 
ues as compiled in Table UTTl 

frequencies in Table ITTTl and included them in Fig. |S|and 
ini Agreement between the experiment and theory is very 
good. 

As a final remark we emphasize that here only high 
symmetry sites have been considered. A mixed occu- 
pation of different sites or off-symmetry adsorption at 
higher coverage can significantly increase the CO stretch 
frequency. As an example we refer to a recent publication 
on adsorption of CO on the Ni(lll) surface presented by 
Eichler ,56j. 



D. Electronic structure 

The electronic structure provides deep insight into the 
interaction between adsorbate and surface. There are 
many papers which deal with electronic structures and 
their influence on the trends in binding energy or geo- 
metrical structure of the CO molecule on metallic surface 
[Tol | . The importance of the interplay between the geo- 
metric and the electronic structure in the understanding 
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Experiment 


Theory 


surface 


site i-'M-co 


(cm-^) 


uc-o (cm^^) 


e (ML) 


Method Ref. 


cal 
'^M-CO 


cal 


Co 


top 




2012 


0.33 


PAIRS 


[63] 


432 


1989 


Ni 


top 




2044 


0.25 RAIRS,HREELS 


[56] 


430 


2041 




hollow 


400 


1873 


0.25 RAIRS,HREELS 


[56] 


353 


1804 


Cu 


top 


346 


2072 




PAIRS 


[64] 


323 


2038 






331 


2077 


0.33 


EELS, PAIRS 


[65] 


- 


- 








2075 low(0.5L) 


PAIRS 


[66] 


- 


- 


Ru 


top 


445 


1980-2080 


0.07 


EELS 


[67] 


439 


1990 






447 


1990 


0.33 


PAIRS 


[68] 


- 


- 








1992 


0.25 


HREELS 


[69] 


- 


- 


Rh 


top 




2000 


0.33 


HREELS 


M 


468 


2029 






480 


1990 


low 


HREELS 


m 


- 


- 








2049 


0.2 


HREELS 


m 


- 


- 




hollow 




1830 


0.2 


HREELS 




328 


1782 


Pd hollow 




1848 


0.33 


PAIRS 


m 


319 


1810 








1823 


<0.25 


PAIRS 


m 


- 


- 








1823-1850 


<0.33 


HREELS 


[73] 


- 


- 




bridge 




1920 


0.5 


PAIRS 


[60] 


348 


1883 


Ag 


top 




2137 


low(2L) 


HREELS 


[74] 


203 


2050 


Ir 


top 




2063 


low(lL) 


PAIRS 


[48j 


505 


2041 








2028-2090 


0-0.71ML 


PAIRS 


[Z4] 


- 


- 








2065 


0.25 


PAIRS 


[76] 


- 


- 


Pt 


top 


470 


2100 


0.24 


EELS 


m 


487 


2081 








2090 




SFG 


[77j 


- 


- 






464 




0.5 


PAIRS 


m 


- 


- 






467 


2104 


0.5 


PAIRS 


[79] 


- 


- 








2093 


0.07 


PAIRS 


[80] 


- 


- 








2095 


0.1 


PAIRS 


[81] 


- 


- 








2100 


0.25 


PAIRS 


[82] 


- 


- 




bridge 


380 


1850 


0.24 


EELS 


L50] 


372 


1880 






376 




0.5 


PAIRS 


[78] 












1855 


0.5 


PAIRS 


[79j 












1858 


0.07 


PAIRS 


m 








hollow 




1810 


0.5 


PAIRS 


[82] 


352 


1793 








1736 


0.51 


PAIRS 


[83] 







TABLE HI: Experimental and calculated metal-CO (z^m-co) and C-0 vibrational frequencies (i^c-o) of CO adsorbed on 3d 
and 4d TM surfaces. The coverage O and references to the experimental studies are added. PAIRS - Reflection- Absorption 
Infra-Red Spectroscopy, (HR)EELS - (High Resolution) Electron Energy Loss Spectra, SFG - Sum Frequency Generation 
spectroscopy. 



of CO adsorption was stressed by Fohlisch et al. In 
our analysis we essentially follow the ideas expressed in 

Refs. j5=m^. 

How is CO adsorbed on the TM surface ? It is gener- 
ally assumed that a major part of the CO - metal inter- 
action can be explained in terms of frontier orbitals (the 
highest occupied molecular orbital (HOMO) and the low- 
est occupied molecular orbital (LUMO) orbitals). The 
Blyholder model is based on the donation from the oc- 
cupied C0-5tT states into empty surface orbitals and the 
back-donation from occupied surface orbitals to the CO- 
27r* orbitals 1]. 

At this point we can divide the literature into three 
groups. To the first belongthose that agree with the Bly- 
holder model completely [H IssL Is^. The second group 
consists of those that agree with Blyholder (in princi- 
ple), but point out that the Blyholder model ignores the 
contribution to the bonding from the 4cr and Itt orbitals. 



Some consider the model not only oversimplified, but pro- 
pose another 2-11* resonance model like Gumhalter et al. 
[901 . And finally to the third group belong those that 
disagree with the Blyholder model and claim that there 
is no back-donation to the 27r* orbital [oil . 

Although in the experimental study of Nilsson et al. 
[92j . the authors emphasize that an atom-specific look 
at the CO adsorption could provide insight into the sur- 
face chemical bond, we have projected the DOS onto the 
molecular CO orbitals. The density of states projected 
onto Itt, 27r*, ha and Aa molecular orbitals (PDOS) of 
a CO molecule in top, bridge, hollow sites and a CO 
molecule far above (~4 A) the surface is depicted in 
the Fig El The weakly bonded CO molecules on the 
Au(lll) surface already render the general interaction 
trends. The corresponding analysis of the orbitals of a 
CO molecule on Pt(lll) can also be found in the paper 
by Kresse et al. [ss! ]. 
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FIG. 10: The projected electronic densities ol states (PDOS) for the CO molecule adsorbed in top, bridge, fee and hep hollow 
sites of the Au (a) and Pt (b) surface. While the upper five panels describe the PDOS for various molecular orbitals of the 
CO molecule, the lower six panels show the PDOS for the substrate atom(s) interacting most with the molecule. The panels 
labelled 'free CO' and 'clean' describe the noninteracting case for comparison (molecule 4 A above the surface). The DOS is 
smoothed by a Gaussian function with a width of 0.2 eV and the Fermi level is located at eV. 



There is no principal difference between the Au and 
Pt DOS. As the CO molecule approaches the surface 
the localized CO orbitals, Scr (not shown) and 4cr, are 
shifted to lower energies, depending on the coordination 
between ~1 and 3 eV. Since they are fully occupied they 
play only a minor role in the interaction with the metal, 
characterized as Pauli repulsion. On the other hand, the 
5(T, Itt and 2t:* peaks broaden and dominate the inter- 
action. Again, the broadening and the shift of the peaks 
increase with coordination: the higher the coordination - 
the greater the shift and the broadening. This behavior is 
more visible for the 5a than for the 2Tr* peak. The posi- 
tion of the 5(7 band peak is almost always lower than the 
Itt peak which is shifted by '-^l to 3 eV. The 5a band is 
higher in energy only for on-top adsorption of CO on the 
Ag(lll) surface, which is partially related to the weak 
interaction. A typical shift for the 5a peak is 3 to 4 eV, 
depending on the metal CO bond strength. Finally, also 
the 2tt* band broadens with increasing coordination and 
is therefore partially shifted below the Fermi level (Epermi 
= eV) as shown in FiglTUl 

The detailed interaction picture on different sites de- 



pends on the symmetry of the surface and adsorbate or- 
bitals. For the top site the major orbital interaction (due 
to symmetry) is 5a (CO) - d^^ (metal), whereas for the 
higher coordinated sites the Itt and 2tt* CO molecular or- 
bitals are more important, interacting with the dxz (resp. 
dyz) and in- plane (dx2_y2, dxy) orbitals of the metal 
atoms. The 5a orbital of the CO orbital adsorbed on-top 
hybridizes with the d^^ states of the metallic substrate 
and shifts to lower energies. The dz2 states broaden and 
split into a 5a — (1^2 bonding contribution far below the 
Fermi-level 7.5 eV) and 5a — d^i anti-bonding contri- 
butions located above the bonding peak and partly even 
above the Fermi-level. This interaction would be repul- 
sive (Pauli like) if the 5a — (1^2 were not pushed partly 
above the Fermi-level. A 5a depletion (donation from CO 
to metal) is in accordance with the Blyholder argument 
of donation of electrons from the adsorbate to the sur- 
face 1]. The interaction between the 5a orbitals and the 
metal s-band is attractive, but depends on the amount of 
electrons accumulated in the newly created molecular or- 
bitals of the CO-metal system. If we consider a depletion 
of the 5a orbital, then the interaction with the metallic 



s-band becomes more attractive. 

The interaction between the l7r(resp. 27t*) orbitals and 
the substrate is more complex jsj; we can distinguish 
four contributions: (i) The main part of the Itt CO or- 
bital is located at energies around ~ — 6 eV. (ii) During 
the adsorption the peak broadens at the higher energy 
end up to the Fermi level. This broad state is usually 
called dfi- orbital. All the interactions with the Itt orbital 
are located below the Fermi level and have repulsive char- 
acter |93|. (iii) Furthermore, in the same energy region 
is a contribution from a partially occupied 27r* orbital 
which increases with coordination and decreases with the 
d-band filling (see Fig. I10|l . This contribution develops 
into a broad peak for Pt. (iv) Finally, at higher energies 
{^3 eV) the anti-bonding 27r'^-dyz(dxz) hybridized orbital 
is visible. 

The later two contributions can be characterized by 
the fractional occupation of the 2tt* orbital, compiled in 
Table Hvl At this point one should remark that our rep- 
resentation of the DOS projected onto molecular orbitals 
is subject to small inaccuracies not in position, but in 
amplitude arriving from the projection of the plane-wave 
components. There exists a tiny hybridization between 
the molecular Itt and 2tt* orbitals, therefore the calcu- 
lated DOS is slightly overestimated at the position of 
the 27T* for the Itt and at the position of the Itt for 2tt*. 
This must be considered when calculating the occupation 
of the 2Tr* orbital. 



Occupation of 2n* 



site 


Co Ni Cu 


Ru Rh Pd Ag 


Ir Pt Au 


top 
bridge 
fee 
hep 


16.0 11.5 8.1 
19.5 14.5 11.1 
29.9 22.1 16.3 
29.7 21.6 16.1 


11.4 10.8 8.8 5.3 

12.5 11.4 10.5 7.0 
18.5 18.5 18.1 11.7 
19.9 19.5 18.2 11.5 


11.4 10.0 6.0 
11.1 10.1 7.4 

20.5 19.8 15.0 
21.4 19.7 14.4 



TABLE IV: Fractional occupation of the 2tt* orbital for the 
CO molecule in top, bridge, fee and hep sites on the close- 
packed transition metal surfaces in %. 



The occupation of the 2Tr* orbital increases with coor- 
dination, but decreases with d-band filling and also from 
3d to the 4d and 5d elements. We can further correlate 
the site preference between the two competitive (fee and 
hep) hollow sites to the occupation of 27r* orbital. As 
the occupation of the orbital for the fee site compared to 
hep site increases, the site preference changes from the 
hep to the fee hollow sites. This effect is probably due to 
a better interaction with the metal atom in the second 
layer below the hep site for metals with smaller lattice 
constants. 



E. Charge density redistribution 

In this section, we illustrate the ideas derived in the 
previous section from the analysis of the DOS by analyz- 
ing the charge density and the difference in the charge 
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(a) 4a (b) 5g (c) 1 71 




(d) 4a (e) 5a (f) 1 n 




FIG. 11: Charge density of (a) 4a, (b) 5a, (e) Itt CO orbitals 
for the free molecule and the difference (Ap = Pco+Au{iii) — 
Pau{iii) — Pco) after the adsorption in the top site of the 
Au(lll) surface (d, e, f). Figures (g) and (h) show the change 
in the metal charge density due to CO adsorption in the en- 
ergy interval (—5.9, —1.0) and just below Fermi level (—0.2, 
Ef). The total charge density difference due to adsorption is 
shown in Fig. (i). Dark regions: charge accumulation, light 
regions: charge depletion. 



density for the molecular orbitals of the CO molecule ad- 
sorbed on-top on Au(lll) surface as shown in Fig. ITTI 
At first we show the Aa, ba and Itt orbital-decomposed 
charge density of the free CO molecule (panels a-c) . We 
can see the accumulation of the charge density for the 4cr 
and Itt orbitals around the oxygen atom and for the 3cr 
(not shown) and 5a orbitals on the carbon atom. Fig. 
IllH . Illb and lilt illustrate what happens to these or- 
bitals after the CO molecule is chemisorbed on-top on 
the Au(lll) surface. Within the 4<t orbital charge den- 
sity moves from the C-0 bond and from above the oxy- 
gen atom to the region below the carbon (Fig. IlllF). 
Fig. Illb shows redistribution of charge in the 5(j molec- 
ular orbital: charge is transferred from the side located 
closer to the metal to the opposite side of the carbon 
and oxygen atom. This charge transfer is in line with the 
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experimental observation of a lo ne p air localized on the 
oxygen atom as reported in Ref . |93| • There is only little 
transfer of charge density from the oxygen orbitals with 
Px/Py character into the metal-carbon bond (Itt). This 
re-distribution could be related to the slight inaccuracy 
of the orthogonalization of Itt and 2n* orbitals. 

The newly formed hybrid states, either by broadening 
of the Itt orbital (d#) or by occupation of 2-k* orbitals is 
shown in panels g and h. The ds- orbital shows a charge 
transfer from the d-orbitals towards the 7r-like molecular 
orbitals. We can even spot a small depletion of charge 
at the metal atom which is only a 2nd nearest neighbor 
of the CO, charge which is attracted towards the metal 
atom interacting with the CO molecule. Furthermore, 
27f* states become occupied as can be seen in Fig. Illh . In 
the panel i, we present the total charge density difference 
of all states, i.e. the sum of all mentioned charge density 
redistributions. 

Summarizing, we can say that there is a depletion of 
charge in the substrate out-of-plane d-orbitals and accu- 
mulation of charge density in the in-plane d-orbitals of 
the substrate. The charge density on CO is re-distributed 
and additionally some charge is transferred into 27r*-like 
orbitals, which weakens the C-0 bond. The establish- 
ment of the M-CO bond is reflected by the accumulation 
of charge between the bonding species. 



V. IMPROVING SITE PREFERENCE 

As previously mentioned the calculations for Cu, Rh, 
Pt surfaces predict a wrong site preference for CO ad- 
sorption and the adsorption energies for different adsorp- 
tion sites on Co, Ag and Au are almost degenerate, so 
that small changes in the setup (k-points and neglect 
of surface relaxation) can alter the site preference. Can 
harder pseudopotentials (lower rcut in the pseudopoten- 
tial generation) or different GGA approaches influence 
the site preference ? In order to shed light onto this 
question we repeated our calculations performed with the 
PW91 exchange-correlation functional and cut-off energy 
Ecut = 450 eV with a harder pseudopotential which re- 
quires a higher Ecut — 700 eV and used two different 
GGA functionals: PW91, RPBE. The use of RPBE is 
known to result in lower adsorption energies for CO, im- 
proving the agreement with experiment in most cases 
at almost unchanged geometrical structure. The results 
for RPBE with Ecut = 700 eV (calculated at the opti- 
mized PW91 geometry) are presented in Fig. 1121 

Switching to the RPBE GGA changes the energetics 
drastically by ^--^0.3 eV and the values agree well with 
the experiments for strong chemisorption. On the other 
hand, the deviation from the experimentally determined 
adsorption energy increases for weak bonding. For Ag 
and Au, the RPBE even predicts CO adsorption to be en- 
dothcrmic. Furthermore, the RPBE functional removes 
the degeneracy of the adsorption sites for Co(OOOl) and 
Ag(lll) surface and on-top sites become clearly favoured. 
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FIG. 12: CO adsorption energies on TM surfaces for the 
RPBE exchange-correlation functional and energy cut-ofF 
700 eV. Experimental CO heat of desorption are labelled by 
the stars (see Table HJ, calculated adsorption energies for top 
sites are presented with triangles, bridge - squares and fee 
hollow - diamonds and hep sites with circles. 



The favoured sites for CO on Rh, Pt, and Cu still con- 
tradict the experimental findings. However, in all these 
cases the corrugation of the potential energy surface is 
very small compared to metals where CO satisfies the 
experimental site preference (e.g. Ir and Pd 3]). This 
also indicates that it is harder for the CO molecule to 
find the optimal adsorption site on these metals. 

To make sure that this large change in Eads is not the 
effect of the Ecut (larger basis set), but of the different 
functional, we have recalculated all systems with a harder 
PW91 pseudopotential which also requires an energy cut- 
off of Ecut = 700 eV. Besides a small reduction of the ad- 
sorption energies (weaker bonding) by about 50 meV the 
results are identical to those obtained for Ecut = 450 eV. 

Summarizing, harder pseudopotentials for C and O 
lead to a slight decrease of the adsorption energies. The 
RPBE exchange-correlation functional reduces the ad- 
sorption energies by ~0.3 eV and improves the differ- 
ences in adsorption energies for different adsorption sites. 
However, the site preference remains the same as for the 
PW91 functional. 



VI. DISCUSSION 

In this section, we discuss trends which can be deduced 
from our results for CO adsorption on closed-packed tran- 
sition and noble metal surfaces. 
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A. Correlations between adsorbate geometry, 
frequency and adsorption energy 

A strong correlation exists between the C-0 bond 
length and the stretching frequency of the CO molecule 
j/C-o- During adsorption the 27r* orbital is partially pop- 
ulated and the C-0 bond is weakened. This filling in- 
creases from the right to the left in the Periodic Table and 
with increasing coordination of the adsorption site and 
hence the CO bond length increases. A weaker and there- 
fore longer C-0 bond implies stronger CO-metal bonding 
and results in lower C-0 stretching frequencies. The de- 
pendence of dc-o on the coordination is smallest for the 
3d and largest for 4d metals as can be seen in Fig. [7| 
The trend goes as follows: the higher the reactivity of 
the surface, the stronger the bonding and population of 
27r "^-derived orbitals, resulting in a more elongated C-0 
bond length and a lower stretching frequency. 
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FIG. 13: Correlation between the calculated CO bond length 
(dc-o) and the C-0 stretching frequency (i/c-o) for all high- 
symmetry adsorption sites. 

In Fig. ^|the CO stretching frequency is plotted ver- 
sus the molecular bond length exhibiting a linear rela- 
tionship between these two properties, which is almost 
independent of the coordination of the CO molecule. The 
linear regression fit is characterized by a correlation co- 
efficient of R = 0.97. 

As expected, pronounced stretching of the CO-bond 
(reflecting a stronger adsorbate-substrate interaction) is 
accompanied by a decrease of the adsorption height. In 
Fig. 1141 we present this linear relationship between the 
vibrational frequency i/qq and the surface-carbon dis- 
tance. Remarkably, a trend established by the noble 
metals alone would be characterized by a much bigger 
correlation coefficient of R = 0.98 (fit is not shown) com- 
pared to the correlation coefficient of R = 0.93 for the 
complete dataset. 

Although we cannot directly relate the height of the 
molecule on the surface to the adsorption energy it serves 



2100 r 



2000 - 



1900 - 



1800 - 



170QL 



• 


topTM 


o 


top noble 


■ 


bridge TM 


□ 


bridge noble 


♦ 


fccTM 


o 


fee noble 


▲ 


hcpTM 


A 


hep noble 




linear fit 



R = 0.93 



1.2 



1.4 



1.6 



1.8 
.[A] 



2.0 



2.2 



FIG. 14: Correlation between the calculated S-C distance 
(ds-c) and the C-O stretching frequency (uc-o) for all high- 
symmetry adsorption sites. 



as a rough estimate of the metal-adsorbate interaction. 
As a consequence of these two linear relationships (dc-o 
to i/QO and ds-c to dc-o) a similar correlation exists 
between the ds-c and dc-o (Fig. IT^ . Again, we have 
divided our dataset into two groups: noble and transition 
metals and we provide separately a linear fit for each 
group of out set. Again, the correlation coefficient is 
higher for noble (Rnobic = 0.98) than for transition metals 
(Rtm = 0.93). 

Above mentioned linear relationships have an impor- 
tant consequence: if we know one of four measured prop- 
erties {vc-o, dc-o, dg-c or d-M-c) we can easily esti- 
mate all other linearly related properties. 
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FIG. 15: Correlation between the calculated CO bond length 
(dc-o) and the S-C distance (ds-c) for all high-symmetry 
on-surface sites. 
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B. Correlation between the electronic properties of 
the substrate and adsorption energy 

As already mentioned above, one of the most widely 
discussed models for molecular adsorption on transition 
metal surfaces is the d-band model of Hammer and 
N0rskov 3, ^ relating the chemisorption energy to the 
position of the metallic d-band relative to the molecular 
orbitals of the adsorbate and to the coupling matrix el- 
ements and overlap integrals between the metal d-states 
and adsorbate states. In a simplified form of the argu- 
ment, the adsorption energy is related to the d-band posi- 
tion only. This simplified d-band model has recently been 
examined by Lu et al. 94] for the adsorption of CO, H2 
and ethylene. This analysis led the authors to a critique 
of the d-band model and of DFT calculations of adsorp- 
tion energies in general. For CO it was pointed out that 
the reduction of the adsorption energies with increasing 
distance of the d-band center from the Fermi energy fol- 
lows qualitatively the predicted trend, but quantitatively 
the theoretically predicted dependence is much stronger 
than what is observed experimentally. 

In Fig. El we present calculated CO adsorption en- 
ergies (PW91) with respect to the d-band centers taken 
from Fig. ^ We can clearly establish a linear relation- 
ship between position of the d-band center and the ad- 
sorption energy for CO adsorption on noble metals. On 
the other hand, the values for the transition metals show 
a large scatter, hence what really determines the trend 
are the noble metals. It does not mean that there is no 
relation between the adsorption energy and the d-band 
center, but it must be corrected for the variation of the 
coupling matrix elements between different metals, as in 
the original argument of Hammer and N0rskov |^ . 
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FIG. 16: Correlation between the calculated position of the 
d-band center (where at least 9.5 electrons are accumulated, 
see Fig. |^ for all high-symmetry on-surface sites. Empty 
symbols are for noble metals and full symbols are for transi- 
tion metals. 



C. Site preference - limitation of DFT 

As a by-product of this study the performance of the 
PW91 functional with different Ecut (harder potential) 
and the RPBE functional (well suited for CO and NO 
adsorption) could be tested on a large number of systems. 

A calculation with harder pseudopotential for the same 
exchange-correlation functional (PW91) leads to small 
changes in the heat of adsorption which is closer to the 
experimental value for TMs, but further for noble metals. 
The site preference does not depend on the basis set for 
PW91 exchange-correlation functional. 

The most important difference between the PW91 and 
RPBE functional with more accurate pseudopotential 
(Ecut — 700 eV) is in the absolute value of the of the 
adsorption energy Ai?ads ^0.3 eV. Moreover, the differ- 
ence in the adsorption energy between previously almost 
degenerate top and hollow adsorption sites is increased 
(Co, Ag, Au). However, all these changes do not correct 
the prediction of a wrong adsorption site for Cu, Rh and 
Pt. 

Hence the question "Why does DFT fail in predict- 
ing the right adsorption site as well as energy ?" arises. 
In the past, it has been proposed that a relativistic 
treatment of transition metals could solve the problem 
[95! 113 ■ However, the results presented to sup- 
port this suggestion are not really convincing, being 
based on rather thin slabs. Also we note that the work- 
function of the transition metals is well described and the 
electronic structure calculations are in good agreement 
with spectroscopic results. A more promising approach 
starts from the observation that 27r*-backdonation is the 
most important contribution to the adsorbate-substrate 
bond. It is well known that DFT calculations tend 
to underestimate the HOMO-LUMO gap in CO, fa- 
cilitating back-donation. Very recently, Kresse et al. 
[ssj l have used an LDA-I-U approach, leading to an in- 
creased HOMO-LUMO splitting due to an additional 
on-site Coulomb repulsion. It was shown that with an 
increasing HOMO-LUMO gap on Pt(lll) 27r* backdo- 
nation favouring hollow-adsorption is reduced and site- 
preference switches to the correct top-site. 

Another way to improve the HOMO-LUMO gap is 
to use hybrid functionals mixing DFT-exchange with 
exact (Hartree-Fock) exchange: for gas-phase CO the 
PW91 functional predicts a HOMO-LUMO gap of 6.8 
eV, whereas the B3LYP hybrid functionals predicts a 
splitting of 10.1 eV. However, B3LYP calculations can 
at the moment only be performed using localized basis 
sets and small clusters. Very recently Gil et al. 
have performed a series of cluster-calculations for CO on 
Pt(lll) using both PW91 and B3LYP functionals and 
slab calculations using PW91. The conclusion is that 
the B3LYP functional reduces the overestimation of the 
stability of the hollow site compared to PW91 for a given 
cluster size. However, the results also show a pronounced 
cluster-size dependence: for small clusters the preference 
for the hollow-site is more pronounced than for larger 



14 



clusters (the largest cluster considered in this study is 
Ptig). For a comparison of cluster and slab calculations 
it is concluded that a slab calculation, if performed with a 
B3LYP functional should favour top adsorption, in agree- 
ment with experiment. 

Current work is devoted to the extension of the 
LDA+U approach to other metals and to the investi- 
gation of the influence of meta-GGA functionals (incor- 
porating a dependence on the kinetic energy density). 

VII. CONCLUSIONS 

The adsorption behavior of the CO molecule on the 

surfaces of closed-packed late transition metals, namely 
Co, Ni, Cu, Ru, Rh, Pd, Pd, Ag, Ir, Pt, Au has been 
analyzed using ab-initio DFT methods. To rationalize 
the trends in the CO adsorption we performed investi- 
gations of the geometric structure (metal surface-C and 
C-0 bond length), the electronic structure (projected 
DOS, charge flow analysis) and the vibrational proper- 
ties (fc-o and i/m-co)- In addition, we examined the 
influence of different GGAs (PW91, RPBE) on the CO 
site preference and adsorption energies. 

We demonstrate that current DFT slightly overesti- 
mates substrate lattice constants (1 to 2%) and under- 
estimates the work-function (by about 0.3 cV). We com- 
pared our calculated values for the geometric structure of 
the adsorbate-substrate complex with experimental data 
and we find them in very nice agreement. Furthermore, 
the C-0 and M-CO stretching frequencies agree well with 



experiments. We illustrate the detailed interaction pic- 
ture of the CO molecules for different adsorption sites on 
TM surfaces by means of density of states projected onto 
molecular CO orbitals and change of the charge density 
due to adsorption. 

Present density functionals overestimate the adsorp- 
tion energy for sites with high metal coordination com- 
pared to sites with low metal coordination. Moreover, the 
underestimated energy difference between HOMO {5a) 
and LUMO orbitals contributes to the CO over- 

binding. DFT fails in predicting the site preference for 
CO on the Cu (by > 30 meV), Rh (by > 30 meV), and Pt 
(by 90 meV) close-packed metal surfaces at 0.25 ML 
coverage. 

Different GGA functionals reduce these discrepancies, 

but neither harder potentials for CO, nor different GGA 
functionals correct the results for CO adsorption on 
these substrates. One should also point out that the 
adsorption energies calculated using RPBE exchange- 
correlation functional give results which arc much closer 
to the experimental values for TM surfaces, but predict 
CO adsorption to be endothermic on the noble metals 
(Ag and Au). 

Recent results have shown that a correct site- 
preference can be achieved by correcting the HOMO- 
LUMO gap using techniques beyond DFT. This is af- 
ter all not entirely surprising, since the analysis demon- 
strates the importance of the 27r*-backdonation effects, 
i.e. electron transfer to excited states of the free 
molecule. 
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